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Summary. Long range charge transfer experiments in DNA oligomers and the 
subsequently measured — and very diverse — transport response of DNA wires in 
solid state experiments exemplifies the need for a thorough theoretical understanding 
of charge migration in DNA-based natural and artificial materials. Here we present a 
review of tight-binding models for DNA conduction which have the intrinsic merit of 
containing more structural information than plain rate-equation models while still 
retaining sufficient detail of the electronic properties. This allows for simulations 
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of transport properties to be more manageable with respect to density functional 
theory methods or correlated first principle algorithms. 

1 Introduction and motivation 

Within the class of biopolymers, DNA is expected to play an outstanding role 
in molecular electronics [1]. This is mainly due to its unique self-assembling 
and self-recognition properties which are essential for its performance as car- 
rier of the genetic code. It is the long-standing hope of many scientists that 
these properties might be further exploited in the design of electronic cir- 
cuits [2-6] . In the last decade of the 20th century transfer experiments in nat- 
ural DNA in solution showed unexpected high charge transfer rates [3,7-10]. 
This suggested that DNA might support charge transport. In contradistinc- 
tion, electrical transport experiments carried out on single DNA molecules 
displayed a variety of possible behaviors: insulating [11, 12], semiconduct- 
ing [13, 14] and ohmic-like [15-18]. This variation might be traced to the high 
sensitivity of charge propagation in DNA to extrinsic (interaction with hard 
substrates, metal-molecule contacts, aqueous environment) as well as intrin- 
sic (dynamical structure fluctuations, base-pair sequence) influences. Recently, 
experiments on single poly(GC) oligomers in aqueous solution [17] as well as 
on single suspended DNA with a more complex base sequence [14] have shown 
unexpectedly high currents of the order of 100-200 nA. Again these results, 
if further confirmed, suggest that DNA molecules may support rather high 
electrical currents given the right environmental condition. 

The theoretical interpretation of these recent experiments and, in a more 
general context, the elucidation of possible mechanisms for charge transport 
in DNA have not, however, been unequivocally successful so far. While ah 
initio calculations [19-28] can give at least in principle a detailed account of 
the electronic and structural properties of DNA, the huge complexity of the 
molecule and the diversity of interactions present preclude a complete treat- 
ment for realistic molecule lengths. When interactions with counter ions and 
hydration shells or vibrational degrees of freedom are to be considered, the 
situation easily becomes untractable. On the other hand, model-based Hamil- 
tonian approaches to DNA [29-44] have been already been discussed in great 
detail and can play a complementary role by addressing single factors that 
influence charge transport in DNA. However, here it is of course clear neither 
a-priori nor a-posteriori (given the aforementioned experimental situation) 
which model should be used. Somewhat mirroring the experimental situation, 
a large variety of models exist and the results are not necessarily consistent 
across different models. 

In this contribution, we review tight-binding models of DNA which have 
been proposed in the literature and argued to reproduce experimental [29] 
as well as ab-initio results [45]. We first concentrate on simple one- and two- 
channel models of DNA in which the main transport mechanism is concen- 
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trated in the effects of 7r-overlap in the base or base-pair sequences. The 
models are usually constructed to take into account the HOMO-LUMO gap 
of the single base pairs similar to many of the DFT-based studies. A main 
feature of the next class of models is the presence of sites which represent 
the sugar-phosphate backbone of DNA but along which no electron transport 
is permissible. These models construct a gap due to transversal perturbation 
of the 7r-stack, i.e. even when the onsite energies are constant. The aim of 
this review is thus to explain the present state of affairs in the tight-binding- 
model-based approach and we will be very brief on many others aspects of 
the charge migration problem, as these are already well treated in the other 
chapters of this book. 

2 The electronic structure of DNA 

DNA is a macro-molecule consisting of repeated stacks of bases formed by 
either AT (TA) or GC (CG) pairs coupled via hydrogen bonds and held in 
the double- helix structure by a sugar-phosphate backbone. In Fig. 1, we show 
a schematic drawing. The electronic energetics of a double-stranded DNA 
chain should take into account three different contributions coming from (i) 
the nucleobase system, (ii) the backbone system and (iii) the environment, 
as it is sketched in Fig. 2. Attending to the energies involved in the differ- 
ent interactions, the resulting energy network can be hierarchically arranged, 
starting from high energy values related to the on-site energies of the bases 
and sugar-phosphate groups (8 — 12 cV) [46,47] passing through intermediate 
energy values related to the hydrogen bonding between Watson-Crick pairs 
(~ 0.5 eV) [46] and the coupling between the bases and the sugar moiety 
(~ 1 eV) [47] and ending up with the aromatic base stacking low energies 
(0.01 — 0.4 eV) [46,48]. The energy scale of environmental effects (1 — 5 eV) is 
related to the presence of counter ions and water molecules, interacting with 
the nucleobases and the backbone by means of hydration, solvation and charge 
transfer processes. It is about one order of magnitude larger than the coupling 
between the complementary bases, and about two orders of magnitude larger 
than the base stacking energies. 

We emphasize that in many of the models to be reviewed later in this chap- 
ter, simplified assumptions about these energy scales are employed. Mostly, 
however, the ionization energies £q = 7.75 eV, ec = 8.87 eV, £a = 8.24 eV 
and €t = 9.14 eV, [48-52] are taken as suitable approximations to the onsite 
energetics at each base. 
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3' end 5 ' end 

Fig. 1. The chemical composition of DNA with the four bases Adenine (A), 
Thymine (T), Cytosine (C), Guanine (G) and the backbone. The backbone is made 
of phosphorylated sugars shown in light grey, the nucleobases are indicated in dark 
grey. 



3 Numerical techniques for charge transport in the 
quantum regime 

Before we in the following turn our attention to the variety of simplified models 
which have been proposed to capture the essential charge transport features 
of DNA, let us briefly recall some of the techniques used to investigate these. 

There are several approaches suitable for studying the transport properties 
of quasi-one-dimensional tight-binding models for long DNA and these can 
be found in the literature on transport in solid state systems, or, perhaps 
more appropriately, quantum wires [53]. Since the variation in the sequence 
of base pairs precludes a general solution, one normally uses methods well- 
known from the theory of disordered systems [54,55]. The main advantage 
of these methods is that they work reliably (i) for the relatively short DNA 
strands ranging from 13 base pairs (as in DFT studies [56]) up to 30 base 
pairs length which are being used in the nanoscopic transport measurements 
[13] as well as (ii) for somewhat longer DNA sequences as modeled in the 
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Fig. 2. Sketch illustrating the overall energetics of a double-stranded DNA chain. 



electron transfer results and (iii) even for complete DNA sequences which 
contain, e.g. for human chromosomes up to 245 million base pairs [57]. We 
measure the effectiveness of the electronic transport by various measures such 
as the localisation length £, participation numbers, etc. These roughly speaking 
parameterise whether an electron is confined to a certain region of the DNA 
(resulting in insulating behavior) or can proceed across the full length L of 
the DNA molecule (metallic behavior). 

3.1 Recursive Green function technique 

The first method one can use is the recursive Green function approach pi- 
oneered by MacKinnon [58,59]. It can be used to calculate the dc and ac 
conductivity tensors and the density of states (DOS) of a d-dimensional dis- 
ordered system and has been adopted to calculate all kinetic linear-transport 
coefficients such as thermoelectric power, thermal conductivity, Peltier coeffi- 
cient and Lorenz number [60,61]. Briefly, the approach utilizes the advanced 
and retarded Green's functions, G~{E — i0 + ) and Q + (E + i0 + ), respectively, 
and the usual definition [(E ± i0 + )<% — Hy] = Sij, where Gfj is the ma- 
trix element (gI^Ij) and Hij is similarly the matrix element of the Hamilto- 
nian [61]. 8^ denotes the Kronecker S between basis states {\i}}. If contains 
only nearest-neighbor connections, then these expressions can be written in- 
teractively as 
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-H ii+1 Gf +lj = S iS - [{E ± »0)<J tf - G% + Hu_iGtu ■ (1) 

Here Hu±i are the terms in the Hamiltonian connecting slice i with its neigh- 
boring slices i ± 1. If we now reinterpret the left index i as a pseudo time, 
then we see that the future Green function slice i + 1 can be constructed 
by the present slice at i and the previous slice at i — 1. The method is well 
suited to study coherent transport properties and can be extended to include 
incoherent processes as well [62]. 

3.2 Transfer and transmission matrix approach 

The next method of choice is the iterative transfer-matrix method (TMM) 
[54, 63-66] which allows us in principle to determine the localisation length 
£ of electronic states in systems with varying cross section M and length 
L > M. This localization length describes the decay of the wave function for 
transport along a quasi one-dimensional system and £ may be used as a rough 
guide of the extend of electronic states. 

For disordered systems, typically a few million sites L 3> M are needed 
to achieve reasonable accuracy for £ [54]. However, in the present situa- 
tion we are interested in finding £ also for DNA strands of typically only 
a few hundred or a few ten thousand base-pair long sequences. Thus in or- 
der to restore the required precision, one modifies the conventional TMM 
and now performs the TMM on a system of fixed length L . This modifi- 
cation has been previously used [67-69] and may be summarized as follows: 
After the usual forward calculation with a global transfer matrix T Lo , we 
add a backward calculation with transfer matrix T r b . This forward-backward- 
multiplication procedure is repeated K times. The effective total number 
of TMM multiplications is L = 2KLo and the global transfer-matrix is 
tl = (7j 7i ) • It can be diagonals as for the standard TMM with K — > oo 

to give t|tl — > exp[diag(4KXo/£*)]- The largest for alii = 1, . . . , M then 
corresponds to the localisation lengths of the electron on the DNA strand 
and will be measured in units of the DNA base-pair spacing (0.34 nm). Let 
us emphasize that the above approach converges even for L < £. However, 
in that case, values of £ clearly are dominated by finite-size and boundary 
effects and their significance is no longer quantitative, but qualitatively in- 
dicates extended states smeared out over the finite system length L. Last, 
the transmission coefficient Tl(E), related to the Landauer conductance g 
via g = (2e 2 /h)T L (E)/(l - T L (E)) [70-72], is defined in terms of the matrix 
elements of tl ■ 

3.3 Attaching leads 

Let us assume that, as a first approximation, we can consider a DNA model 
in terms of a linear chain with an orbital per site, where each lattice site 
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represents a base pair. The ends of the chain are connected to leads mod- 
eled as semi-infinite one-dimensional chains of atoms with one orbital per 
site. Broadly speaking, one expects the binding to metallic leads would af- 
fect the electronic structure of the molecule. If so, we should consider the 
states belonging to the coupled molecular-metallic system rather than those 
of the molecular subsystem alone [73] . Thus we shall consider henceforth that 
the coupling between the contacts and the molecule is weak enough, so that 
the lead-molecule-lead junction can be properly described in terms of three 
non- interacting subsystems [74,75], according to 

H = # D NA - Contact (|0)<1| + \N)(N + 1| + h.C.) 

— oo 

+ ^£Lcad|fc)(fc| -*Lcad|fc- + h.C. 

k=0 
+oo 

+ Yl £ Lead|0<*l -t Lead \l)(l + l\ + h.C. (2) 
k=N+l 

In Eq.(2), #dna is the DNA Hamiltonian, the second term describes the 
DNA-lead contact, and the last two terms describe the contacts at both sides 
of the DNA chain, where N is the number of base pairs, e n are the on-site 
energies of the base pairs, icontact is the hopping strength between the leads 
and the end nucleotides, £Load is the leads on-site energy and tLcad is their 
hopping term. 

The Green function methods are well-suited to include contact effects since 
their boundary conditions at the contacts require specification of a suitable 
Green function in the leads which can be chosen to model the geometry of 
contacts. The TMM usually starts assuming a particle- like injection of carriers 
into the transport channels and a proper treatments of leads is lacking, but the 
extracted localization lengths at least for long chains are largely independent 
of the exact choice. Irrespective of numerical methods used, most earlier tight- 
binding studies assumed perfect coupling to metallic leads or simply ignore 
the issue altogether. The role of contact effects within the TMM framework 
was recently reported for a poly(GACT) tetra-nucleotide in Ref. [76] in terms 
of two contact matrices which explicitly take into account the presence of the 
^Contact hopping integral. Depending on the value adopted for ^contact , the ob- 
tained transmission coefficient does not reach in general the full transmission 
condition Tl(E) = 1 due to the symmetry breaking related to the coupling 
of the G (T) end nucleotides at the left (right) leads, respectively. This ex- 
treme sensitivity is due to interference effects between the DNA energy levels 
and the electronic structure of the leads at the metal-DNA interface, and in- 
dicates that the optimal system configuration for efficient charge transfer is 
determined by the resonance condition icontact = y/i ■ iLcad- Quite interest- 
ingly, one realizes that, due to resonance effects, a stronger coupling to the 
leads not always result in a larger conductance through the system, in agree- 
ment with the results obtained by Guo and co-workers for the transmission 
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coefficient of poly(G)-poly(C) molecules making use of the Green function 
technique [77] . Subsequent works have exploited the existence of this optimal 
charge injection condition to study the charge migration efficiency through 
more realistic duplex chains (see the contribution by Apalkov, Wang, and 
Chakraborty in this volume) 

In general, modeling the geometry and bonding character of the contact 
at the interface is a very delicate issue, since detailed information on both the 
metal geometry and DNA chemical bonding at the contacts is poorly known 
to date. Consequently, in most modeling of the DNA-contact interface, the 
parameter ^contact deals with the tunneling probability between the frontier 
orbitals, thus roughly encompassing bonding effects at the interface. Recent 
transport experiments have shown that deliberate chemical bonding between 
DNA and electrodes is a prerequisite for achieving reproducible conductivity 
results [12,78,79]. Accordingly, the study of contact effects on the charge 
migration efficiency is an important issue to be considered in realistic models 
of DNA transport. 

4 Tight-binding model approaches 

The ab initio methods are clearly very powerful. However, from a physics 
perspective, the question immediately arises if an even simpler, effective model 
approach, might capture the essentials of charge migration equally well. This 
strategy is known as the tight-binding approach to DNA — note that in 
this language the term tight-binding is used somewhat differently from the 
terminology of theoretical chemistry. This simple approach has been used 
right from the start of the physics involvement in DNA research. The idea is 
to capture the main path-ways of charge migration along the DNA molecule 
stack in a simple model of sites and hopping strengths. Charge transport along 
this model is then described by simple tight-binding orbitals on the sites and 
suitably parametrized hopping onto neighboring sites. The advantage is this 
approach should be clear: once the appropriate onsite energies and hopping 
strengths are known, much larger system sizes can be studied than with the 
ab initio methods. The downside of course is that the determination of the 
effective parameters and in particular the choice which to leave out completely 
will be at least to some degree a matter of personal preference and thus open 
to criticism. 

4.1 Importance of the DNA sequence: one-dimensional models 

The simplest TB model of the DNA stack can be constructed as a one- 
dimensional model as given in Fig. 3. There is a single central conduction 
channel in which individual sites represent a base-pair. Every link between 
sites implies the presence of a hopping amplitude. The Hamiltonian for this 
wire model (i?w) is given by 
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Fig. 3. The line model for electronic transport along DNA corresponding to the 
Hamiltonian given in Eq. (3). Lines denote hopping amplitudes and circles give the 
central (grey) nucleobase pairs. 

L 

H w = ^2-U\i){i + l\-U-i\i){i-M+Ei\i){i\ (3) 

i=l 

where ti is the hopping between nearest-neighbour sites i, i + 1 along the 
central branch and we denote the onsite energy at each site along the central 
branch by e,. L is the number of sites/bases in the sequence. For constant 
ti = t, €i = and open boundary conditions, the spectrum of the model is 
given by 

„__„„(_£.) (4, 

with k = 1, 2, . . . L. For random choice of onsite energies or hopping strengths, 
this model is well-known as the Anderson model [80] with diagonal or off- 
diagonal disorder and its transport properties are governed by one-dimensional 
Anderson localization [81]. 

In order to use this Hamiltonian to model DNA, one needs to know the 
appropriate parameters for onsite energies and hopping strengths [48-52], or, 
alternatively, argues that mostly the statistical properties of these quantities 
determine the transport. For natural DNA sequences, a useful choice for the 
onsite energies might be the ionization potentials mentioned in section 2. But 
since base-pairs are modelled by a single site, the DNA is effectively described 
as a sequence of GC (identical to CG) and AT (or TA) pairs with links between 
like (GC-GC or AT-AT) or unlike (GC-AT, AT-GC) pairs. Thus the model 
parameters for the pairs should be computed as suitable estimates based on 
the ionization potentials of individual bases [48-52] . 

Already such a simple model as (3) allows to study various aspects 
of charge transport in DNA. Electrical transport through invididual DNA 
molecules was studied in Ref. [82], using poly(G)-poly(C) DNA. Individual 
molecules are coupled to external baths [83] , thus leading to partial decoher- 
ence. Good agreement with the experimental results of Ref. [13] is demon- 
strated. A twist angle in the hopping parts of (3) is used in Ref. [84] to model 
the effect of thermal fluctuations on transport in DNA. The participation 
ratio is used to estimate the extent of the electronic states. Assuming that 
inelastic effects due to temperature can be ignored, the paper then computes 
the temperature dependence of the conductivity. The transmission spectrum 
for a chain of poly(G)-poly(C) DNA molecules was studied in Ref. [77] where 




base pairs 
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also disorder and contact effects were taken into account. The model con- 
tains various parameters according to the HOMO/LUMO structure of DNA. 
Furthermore, charging effects, i.e. Coulomb blockade, were studied within a 
mean-field approach. For a DNA chain consisting of AT and GC pairs, Ref. [85] 
investigates structural and dynamical disorder. Here, in addition to the on- 
site energies in (3), also the hopping elements U are chosen according to the 
specific DNA sequence, which itself, however, consists of random sequences of 
A,T,G,C nucleotides. It is shown that both types of disorder can significantly 
influence the transport properties. Ref. [86] studies both (quasi-)coherent and 
incoherent transport regimes using Landauer and Kubo formalism via a con- 
tinued fraction approach for poly(G)-poly(C) and also poly(A)-poly(T) DNA 
chains. Superexchange-like exponential length dependence is found for the 
coherent and Ohmic-likc behavior for the incoherent regimes. 

The next group of studies focuses on the influence that possible correlations 
in both artificial and natural DNA sequences might have on the transport. 
Natural A-phage DNA has been investigated in Ref. [33] within a transfer- 
matrix approach. Transmission spectra are shown to be very different from 
poly(G)-poly(C) DNA. The results are argued to be roughly consistent with 
those from electron transfer studies. The influence of long-range correlations 
in DNA sequences is studied in Ref. [32] . Natural DNA of the first completely 
sequenced human chromosome 22 (Ch22) is compared to articial sequences 
such as random and Fibonacci sequences. Is is found that long-range corre- 
lations induce coherent charge transfer over longer lengths scales, at least for 
Ch22. Ref. [87] uses the same numerical method as Ref. [32] and corroborates 
the results for Ch22 by comparing to a Rudin-Shapiro sequence. An intrigu- 
ing relation between the length of a region in coding DNA versus non-coding 
DNA and a repeatedly higher transport characteristic in coding DNA was 
reported in Refs. [88,89]. 

The influence of temperature and associated structural fluctuations of 
DNA and thus the onsite and hopping parameters has been studied in the 
next group of papers. Ref. [90] investigates a polaronic model in which the 
hopping elements are influenced by vibrations along the chain. It is shown that 
a polaron can indeed form for reasonable values of parameters. Ref. [91] studies 
a similar situation, but also takes into account the rotation between base pairs 
along the DNA stack. The paper is actually aimed at charge transfer and pro- 
poses that the thermal fluctuations are the limiting step for site-to-site charge 
transfer. Polarons, which have a "twist" and can thus model the double-helix 
structure of DNA better, are investigated in Ref. [92]. Non-linear effects are 
taken into account and it is shown that these lead to different polaronic behav- 
ior. In Ref. [93], it is argued that polaronic transport can be trapped by the 
thermal denaturation of poly(G)-poly(C) DNA. Thermal effects are modelled 
by an anharmonic Morse potential. Semi-empirical quantum-chemical calcu- 
lations were performed in Ref. [94] for poly(G)-poly(C) and also for poly(A)- 
poly(T) DNA using a polaron model. Localization lengths of charge states 
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larger than 2000 base pairs have been computed and it was shown that signif- 
icant differences between poly(G)-poly(C) and poly(A)-poly(T) DNA exist. 

Besides temperature, the solution in which DNA is prepared or measured, 
its geometry and bend as well as the properties of the contacts to external 
leads will influence the measured transport characteristics. The influence of 
disorder for (3) has been investigated in Ref. [85]. Contact effects were studied 
by Ref. [76] for poly(GACT) chains. Resonance conditions were identified 
which show that a strong coupling to the leads does not always result in 
larger conductance. 

The simple wire model (3) has also been used for studies of charge transfer. 
Briefly, DNA bridges, containing only AT base pairs, were investigated in 
Ref. [95] and decay lengths comparable with single-step tunneling were found. 
The presence of Kondo bound states [96] leads to long tunneling lengths above 
lOOnm. Similarly, time-dependent random hopping strcnghts were studied in 
Ref. [97] and analyzed in a charge transfer context. Last, a soliton-based 
explanation for charge transfer in long segments of DNA was given in Ref. [98] . 



4.2 Importance of base-pairing: two-channel models 

A central simplification of the wire model is the description of a DNA base- 
pair as a single site. By doing so, one looses the distinction between a pair 
with G (or A) on the 5' end of the DNA and a C (or T) on the 3' side and 
one where C sits on the 5' and G on the 3', i.e. GC is equal to CG. This 
distinction becomes important when considering hopping between base-pairs, 
e.g. the hopping from GC to AT is different from CG to AT because of the 
different size of the DNA bases and thus the different overlap between G to 
A and C to A (and similarly for C to T and G to T) [99]. Furthermore, the 
relevant electronic states of DNA (highest-occupied and lowest-unoccupied 
molecular orbitals with and without an additional electron) are localised on 
one of the bases of a pair only [100]. The reduction of the DNA base-pair 
architecture into a single site per pair, as in the wire model (3), is obviously 
a highly simplified approach. 

This deficiency of the wire model may be overcome by modelling each 
DNA base as an independent site. The hydrogen-bonding between base-pair 
is then described as an additional hopping perpendicular to the DNA stack as 
shown in Fig. 4. There are two central branches, linked with one another, with 
interconnected sites where each represents a complete base. This two-channel 
model is a planar projection of the structure of the DNA with its double-helix 
unwound, and still without regard for the backbone. We note that results for 
electron transfer also suggest that the transfer proceeds preferentially down 
one strand [3] . The Hamiltonian now reads 



(ti,T\i,r)(i + 1,t| +e i>T \i,T){i,T\) 
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T=2 
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Fig. 4. The two-channel model for electronic transport along DNA. The model 
corresponds to the Hamiltonian (5). Electronic pathways are shown as lines, whereas 
the nucleobases are given as (grey) circles. 



where ti yT is the hopping amplitude between sites along each branch r = 1, 2 
and £i jT is the corresponding onsite potential energy. The new parameter t\2 
represents the hopping between the two central branches, i.e., perpendicular 
to the direction of conduction. As before, we may now attempt to use ab- 
initio methods to compute t\i or simply model it relative to the strength of 
the parallel hopping f i]T . For the ordered system with ti iT = t and e^ T = 0, 
the two channel model is just a special case of the 2D rectangular system 
with spectrum — 2t x cos ( j^+l ) — 2t y cos ^ L 7rfc | " 1 ^ , k x = 1,2, ...,L X , k y = 
1,2,..., L y , . Thus we find 



where the minus (plus) sign correspondes to even, -0n,i = i>n,2, (odd, -0n,i = 
— i>n,2) states with the same (opposite) sign for the wave function on each 
strand. For random onsite disorder, the system is again localized and the 
localization lengths are known for different energies and disorder values [101]. 

Iguchi was one of the early authors to suggest that a two-leg ladder model 
might be a useful starting point [102]. A band gap like behavior was found 
in Ref. [103], which also considered the Coulomb repulsion between different 
bases. It was further shown that for engineered DNA — modelled as frustra- 
tion — the band vanishes. Ref. [34] used the two-leg ladder model to study 
the spatial extent of electronic states in long DNA chains. They found that 
the extent varies considerably depending on the sequence, but remains rather 
small. Recently, Caetano and Schulz found very large participation ratios in 
the two-leg ladder at finite system sizes [104]. They speculated that this might 
indicate a transition to effectively delocalized states. But this claim is not ex- 
pected to hold for longer chains [105-107]. The influence of electronic spin and 



+ ii )2 |i,l)(i,2| +h.c. 



(5) 
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interactions has been studied in Ref. [36]. This work concentrates on charge 
transfer aspects and shows that interaction opens a gap in the electronic 
states of AT and GC pairs. Further transport properties of Ch22, as well as 
A-phage and the histone protein, are investigated in Ref. [41] and compared to 
artificial DNA. It is notable that while the model used in [41] is a two-leg lad- 
der, the rungs of the ladder are now modeling not the tt transport channels 
but rather the charge migration along the sugar-phosphate backbone. This 
approach is similar to Ref. [108]. Discrete breather-type solutions caused by 
environmental effects were studied in a two-leg ladder already in Ref. [109]. 
A Morse potential was used to represent hydrogen bonding. The breathers 
were found to be pinned by the discrete lattice or trapped in defect regions. 
A similar model based on the non-linear Schrodingcr equation was studied 
in Ref. [110], where the transport of the solitons was assumed to propagate 
along the sugar-phosphate backbone. 



4.3 Backbone effects: The fishbone model 



This fishbone model, shown in Fig. 5, retains the central conduction chan- 
nel in which individual sites represent a base-pair. However, these are now 
interconnected and further linked to upper and lower sites, representing the 
backbone. The backbone sites themselves are not interconnected along the 
backbone. Every link between sites implies the presence of a hopping ampli- 
tude. The Hamiltonian for the fishbone model (Hp) is given by 

^ = EE {-m{i+i\-ti\i, q )(j,\ 

i=i 9=T4 

+e i \i){i\+e q i \i,q){i,q\) + h.c. (7) 

where ti is the hopping along the central branch and tf with q =f, J, gives the 
hopping from each site on the central branch to the upper and lower backbone 
respectively. We denote the onsite energy at each site along the central branch 
by Si and, additionally, the onsite energy at the sites of the upper and lower 
backbone is given by e|, with q =tl- L is the number of sites/bases in the 
sequence. It is easy to see that the existence of the backbone leads to an 
effectively renormalized and energy-dependent disorder 



at each base pair n on the tt stack. If, we as before for the wire model (3), 
consider the ordered situation ti = t, = t^, et = ef — for a —1, 1, we find 
that the energies are now given by 
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backbone 



base pairs 



backbone 



Fig. 5. The fishbone model for electronic transport along DNA corresponding to 
the Hamiltonian given in Eq. (7). Lines denote hopping amplitudes and circles give 
the central (grey) nucleobase pairs and backbone (open) sites. 



for k = 1, 2, . . . , L. Hence, there is a highly degenerate state at E = corre- 
sponding to all the backbone sites and the original single-band of (4) splits 
into two cosine bands such that 



' - \A 2 + '21 -. I ; \//" + 2/ 



u 



t - \ft 2 + 2tf 2 , t + \Jt 2 + 2t^ 



(10) 

In Ref. [29] it had been shown that this model when applied to an artifi- 
cial sequence of repeated GC base pairs, poly(G)-poly(C) DNA, reproduces 
experimental data current-voltage measurements when tj = 0.37 eV and 
t\ = 0.74 eV are being used. Therefore, we will assume t\ = 2ij and set 
the energy scale by U = 1 for hopping between GC pairs. Furthermore, since 
the energetic differences in the adiabatic electron affinities of the bases are 
small [111], we choose £j = for all i. 

The physics of the fishbone model was first discussed for poly(G)-poly(C) 
in Ref. [29]. In fact, the central sites of the fishbone are to model the G 
nucleotide only, with the effect of the C bases neglected as not so relevant for 
transport due to their different onsitc HOMO/LUMO energies. The model 
was then independently studied by Zhong [112] for random and natural DNA 
sequences and he also found an interesting transport enhancing effect in the 
band gap upon increasing potential disorder. A further study [40] revealed 
that the extent of electronic states in the two bands of the model can be up 
to a few dozen base pairs large. Furthermore, upon adding binary disorder, 
intended to model adhesion of ions from the ionic solution in which DNA 
strands exist, the band gap closes and the size of initially very well localized 
band-gap states can be made to increase substantially [113]. This effect was 
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also studied in Refs. [38,39] where the system was coupled to a phonon bath. 
Here, the band gap was shown to close with increasing temperature and the 
temperature dependence of the charge transmission near the Fermi energy is 
exponential. 

4.4 Backbone in a ladder 

Combining the advantages of the fishbone and the two-channel models, we 
now model each base as a distinct site where the base pair is then weakly 
coupled by the hydrogen bonds. The resulting ladder model is shown in Fig. 6. 
There are two central branches, linked with one another, with interconnected 



x=l 



T=2 




Fig. 6. The ladder model for electronic transport along DNA. The model corre- 
sponds to the Hamiltonian (5) and the reader should compare the figure to Figs. (4) 
and (5). 

sites where each represents a complete base and which are additionally linked 
to the upper and lower backbone sites. The backbone sites as in the fishbone 
model are not interconnected. In fact, first principle calculations, showing that 
the phosphate molecular orbitals are systematically below the base related 
ones, do not favor the possible hopping of charge carriers between successive 
phosphate groups along the backbone. [114] The Hamiltonian for the ladder 
model is given by 
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L 



H L = Y^ ( t hr\i,T){i + l,r| +Si, T \i,T){i,T\) 

i=l Lr=l,2 

+ (t!\i,r)(i,q(r)\+e!\i,q)(i,q\) 
?=T4 



+ *i,2|*,l><*,2| 



+ h.c. 



(11) 



where as before in (5) ti iT is the hopping amplitude between sites along each 
branch r = 1, 2 and is the corresponding onsite potential energy. t\ and 
and ef as in (7) give hopping amplitudes and onsite energies at the backbone 
sites. Also, q(r) =T,| for t = 1,2, respectively. The parameter t\2 represents 
the hopping between the two central branches as for the two channel model 
(5). 

For the ordered system with ti iT = t, = fl, e^ T = e? = 0, we find again 
that the presence of the backbone sites leads to an effective renormalization of 
onsite energies along the two base pair strands with energy-dependent disorder 



t° 



and (t, a) = (1, |) or (2, J.). The energies for even states are 

irk + 



(12) 



E+ = 



1 



-t\ 2 — 2t cos 



L + l 



± 



tl2 + 2£ cos 



irk+ 
L + l 



1 2 



with /c + = 1, 2, . . . , L. Similarly, the odd states have energies 

irk 



(13) 



t\.2 + 2t cos 



L + l 



± 



ti.2 — 2i cos 



irk 
L + l 



(14) 

and fc = 1, 2, . . . , L. Thus we again have two energy bands, with a slightly 
smaller gap, given as 



E G 



U 



t + 



- t + 



t 



1,2 
*1,2 



*1 



*1,S 



2 

*1,2 



*1 



+ \ \ t + 



1-1,2 



3) 



In Ref. [40], electronic transport in this model was measured by the localisa- 
tion length £, which roughly speaking parametrises whether an electron is con- 
fined to a certain region £ of the DNA (insulating behaviour) or can proceed 
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across the full length L (< £) of the DNA molecule (metallic behaviour). Vari- 
ous types of disorder, including random potentials, were employed to account 
for different real environments. Calculations were performed on poly(dG)- 
poly(dC), telomeric-DNA, random- ATGC DNA and A-DNA. The authors 
find that random and A-DNA have localisation lengths allowing for electron 
motion among a few dozen base pairs only. A enhancement of localisation 
lengths similar to the fishbone model (7) was observed at particular energies 
for an increasing binary backbone disorder. In Refs. [100, 115], the model was 
used to study differences in different natural and artificial DNA sequences. 
Specifically, promoter sequences and sequences known to be repetitive from a 
biological point of view were investigated to see whether there were statisti- 
cally relevant differences. Using the same sequences as Ref. [89], no support 
for larger £ values in regions of coding DNA was found. 

5 Conclusions 

In this chapter, we have aimed at giving a review of current models used for 
a simplified, tight-binding-based analysis of charge transport in DNA. While 
the models can be roughly classified according to their geometrical structure, 
many of the presently available results appear somewhat disjoint and are 
nearly as widely spread as in the experimental situation. Let us nevertheless 
attempt to identify some common themes. The vast majority of studies pre- 
sented here agrees that the transport properties upon including some degree 
of energetic disorder — be it strictly random or according to some suitable, 
naturally occurring sequence — tend towards the insulating side. Neverthe- 
less, the size of the electronic states for finite DNA strands might be larger 
and even exceed the distance between contacts. In such a situation, the ex- 
perimental results might find finite currents. This finding seems to be largely 
independent of the set of on site energies and hopping strengths chosen. Also, 
most studies agree that there are differences between natural DNA sequences 
and random DNA with the same ATGC content. However, it is not clear if 
these differences are due to the special choice of DNA strands and simply 
become statistically irrelevant when other DNA sequences are considered as 
well. Thus, a clear correlation between charge transport and a particular DNA 
sequence or parts thereof is yet to be discovered. We emphasize, however, that 
if such a correlation were to be found, we would find it useful if it persists 
across most models reviewed here. 
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